Lattice gauge theory model for graphene 
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The effects of the electromagnetic (e.m.) electron-electron interactions in half-filled graphene 
are investigated in terms of a lattice gauge theory model. By using exact Renormalization Group 
methods and lattice Ward Identities, we show that the e.m. interactions amplify the responses to 
the excitonic pairings associated to a Kekule distortion and to a charge density wave. The effect of 
the electronic repulsion on the Peierls-Kekule instability, usually neglected, is evaluated by deriving 
an exact non-BCS gap equation, from which we find evidence that strong e.m. interactions among 
electrons facilitate the spontaneous distortion of the lattice and the opening of a gap. 
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Graphene, a monocrystalline graphitic film that has 
been recently experimentally realized [1], has highly un- 
usual electronic properties, due to the quasi-relativistic 
nature of its charge carriers [2] . In the presence of long- 
ranged interactions, graphene provides an ideal labora- 
tory for simulating Quantum Field Theory (QFT) mod- 
els at low energies and to possibly observe phenomena 
like spontaneous chiral symmetry breaking and mass gen- 
eration. For a few years, there was essentially no ex- 
perimental signature of electron-electron interactions in 
graphene, due to substrate-induced perturbations that 
obscured their effects. However, the realization of sus- 
pended graphene samples is allowing people to collect 
increasing evidence of interaction effects [3] . On the the- 
oretical side, understanding the properties of a system of 
interacting fermions on the honeycomb lattice is a chal- 
lenging problem, similar to quantum electrodynamics, 
but with some peculiar differences that make its study 
new and non-trivial [4] ; its comprehension is essential for 
graphene and, at the same time, it has relevance for other 
planar condensed matter systems, like high Tc supercon- 
ductors, and even for basic questions in QFT. 

Most theoretical analyses deal with a simplified effec- 
tive continuum model of massless Dirac fermions with 
static Coulomb interactions. Early results based on low- 
est order perturbation theory predicted a growth of the 
effective Fermi velocity v{k.) close to the Fermi points and 
excluded the spontaneous formation of a gap at weak cou- 
pling [5] . The absence of a gap is a serious drawback for 
possible technological applications of graphene. There- 
fore, people started to investigate possible mechanisms 
for its generation. One way to induce a gap is by adding 
an interaction with an external periodic field, which can 
be generated, e.g., by the presence of a substrate [7]. In 
the absence of a substrate, Ref.[8] proposed a mecha- 
nism similar to the one at the basis of spontaneous chiral 
symmetry breaking in strongly-coupled QED3; the appli- 
cability of these proposals to real graphene is a delicate 
issue, due to the uncontrolled approximations related to 
a large-N expansions. Another possible mechanism for 



gap generation is based on a Peierls-Kekule distortion of 
the honeycomb lattice, which is a prerequisite for electron 
fractionalization [9, 10]. One optimizes over the distor- 
tion pattern, by minimizing the corresponding electronic 
energy. In the absence of interactions, a rather strong in- 
teraction with the classical phonon field is needed for the 
formation of a non-trivial distortion [9] , while the effects 
of the electron interactions is still poorly understood. 



In this Communication, we investigate the effects of 
electronic interactions in terms of a lattice gauge the- 
ory model. We show that the electronic interactions am- 
plify the response functions associated to Kckulc (K) or 
Charge Density Wave (CDW) pairings. Moreover, we de- 
rive the exact form of the Peierls-Kekule gap equation in 
the presence of the e.m. electron-electron interactions, 
from which we find evidence that strong e.m. interac- 
tions enhance the Peierls-Kekule instability (despite the 
growth of the Fermi velocity [11], which apparently op- 
poses this effect). The model is analyzed by the methods 
of Constructive QFT, which have already proved effective 
in obtaining rigorous non-perturbative results in many 
similar problems [12]; we rely neither on the effective 
Dirac description nor on a large-N expansion. 



The charge carries in graphene are described by tight 
binding electrons on a honeycomb lattice coupled to a 
three-dimensional (3D) quantum e.m. field. We in- 
troduce creation and annihilation fermionic operators 
VjJ = (ai ,6t F ) = \B\-^ fr^„dki^i e^'^^ for elec- 
trons with spin index a =ti sitting at the sites of 
the two triangular sublattices A^ and As of a honey- 
comb lattice; we assume that A^ = A has basis vec- 
tors li^2 = 5(3,±\/3) and that Ab = Aa + Sj, with 
Si — (1,0) and 62,3 — ^(— 1,±-\/3) the nearest neigh- 
bor vectors; B is the first Brillouin zone. The honey- 
comb lattice is embedded in M. and belongs to the plane 
X3 = 0. The grand-canonical Hamiltonian at half-filling 



is H = Ho + He + Ha, where 
Ho = -ty Yai b- ^ e'-IoSrA(s+sS,,o)ds 

" Z^ /-^ 2;,CT 3+Sj,a 

(1) 
with t the hopping parameter and e the electric charge; 
the couphng with the e.m. field is obtained via the Peierls 
substitution. Moreover, if ng = ^^ ^t cr'^s a (I'^sp- ns = 
J2a ^t,a^x,a) ^r X G A^ (rcsp. X € A^), 
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where ipp is an ultraviolet regularized version of the static 
Coulomb potential. Finally, Ha is the energy (in the 
presence of an ultraviolet cutoff) of the 3D photon field 
A = {A, A^) in the Coulomb gauge. We fix units so that 
the speed of light c — \ and the free Fermi velocity v = 
|i <C 1. If we allow distortions of the honeycomb lattice, 
the hopping becomes a function of the bond length (.^j 
that, for small deformations, can be approximated by 
the linear function t^j = t + 4>x,j, with <l)gj = g{ix,j — i) 
and I the equilibrium length of the bonds. The Kekule 
dimerization pattern corresponds to, see Fig.l, 

(l^xj = 00 + Aq cos {p^{Sj - (5jo - x)) , (2) 

with jo G {li2,3}. In order to investigate the effect 




FIG. 1: The Kekule pattern; the hopping parameter is t + Ao 
and t — Ao/2 on the double and single bonds respectively. 

of the e.m. interactions and the Peierls-Kekule instabil- 
ity, we use the following strategy: we first compute the 
ground state energy and the correlations for (j)gj = by 
exact Renormalization Group (RG) methods [13], find- 
ing, in agreement with previous analyses [6], that the 
quasi-particle weight vanishes at the Fermi points and 
the effective Fermi velocity tends to the speed of light as 
power laws with non-universal critical exponents. In ad- 
dition, the analysis of the response functions and of the 
corresponding exponents indicates a tendency towards 
excitonic pairing; the mass terms of K or CDW type are 
strongly amplified by the interactions. Next, we compute 
the electronic correlations in the presence of a non-trivial 
lattice distortion (pgj and we show that a Kekule dimer- 
ization pattern Eq.(2) is a stationary point of the total 
energy (i.e., the sum of the elastic energy and the elec- 
tronic energy in the Born-Oppenheimer approximation). 



We start with (j>gj = 0; the analysis is very similar 
to the one performed in the continuum Dirac approxi- 
mation in [14] (which we refer to for more details), the 
main difference being that the present lattice gauge the- 
ory model is automatically ultraviolet-finite and gauge 
invariant: this avoids the need for an ultraviolet reg- 
ularization, which can lead to well known ambiguities 
[15]. The n-points imaginary-time correlations can be 
obtained by the generating functional 

gW(j,A)^ /p(dV') /p(dA)e^(^+-^''^)+('^'^) (3) 



where: Vk £^re Grassmann variables, with k = {kQ,k) 
and fco the Matsubara frequency, P{dip) is the fermionic 
gaussian integration with inverse propagator 



-i(k) = - 



ika vQ* (k) 
vfl{k) iko 



(4) 



with Z = 1 and n{k) = ^T,j=i.2,3^'''^^'^^'^ (i^ot"^ that 
^(k) is singular only at the Fermi points k — k^ — 
(0,^,±^)); if ^ = 0,1,2, A^(p) are gaussian vari- 
ables with propagator w^^(p) = Sf,„ J ^ ^^_^^p2^^ , 
where x is an ultraviolet cutoff function; finally V = 
eZ JIjqAq + vjA] + h.o.t., where h.o.t. indicates higher 
order interaction terms in A produced by the Taylor ex- 
pansion of the exponential in Hq and j^ is the bare lattice 
current [18]. We compute W{J,X) via a rigorous Wilso- 
nian RG scheme, writing the fields ip,A as sums of fields 
^(fc)^yi(fe)^ living on momentum scales |k— k^|, |p| ~ Af*^, 
with A; < a scale label and A/ > 1 a scaling parameter; 
the iterative integration of the fields on scales h < k < 
leads to an effective theory similar to (3) with an ultra- 
violet cut-off around the Fermi points of width M^ and 
with effective scale dependent wave function renormaliza- 
tion Zh, Fermi velocity Vh and effective charge Ch- This 
approach works only if e^ does not flow to strong cou- 
pling; the boundedness of eh follows from an exact Ward 
Identity (WI), derived by the lattice phase transforma- 
tion 7/-±, ^ e±'^"-V;±, in 11^^(0, A) [19]: 

p^A(,'')(k,p) = e[S^l^Toik,p) - Toik,p)S^'] (5) 

where ^k is the interacting propagator, To{k,p} = 
and AJ, -^(k, p) is the vertex function [20] 











that, if computed at external momenta |k — kr 



A/'' 



and IpI <^ M^, is proportional to et (if M = 0) o^' ^hVh 
(if ^ = 1,2). Using Eq.(5), we find that en — > e_oo = 
e -I- e^F{e), with F{e) a series in e^ with bounded coef- 
ficients, i.e., the effective charge tends to a line of fixed 
points; F[e) is vanishing at lowest order, see Fig. 2, but 
the WI does not exclude that F(e) is non-zero at higher 
orders. A similar WI implies that the photon remain 
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FIG. 2: The second order graphs contributing to the dressed 
charge e_oo, given by the contribution of the vertex part mi- 
nus the graphs coming from the wave function renormaliza- 
tion (if ^ = 0) or the velocity renormahzation (/i = 1, 2); their 
sum is exactly vanishing, in agreement with the WI Eq.(5). 



massless and, as an outcome of the above procedure, we 
get an expansion of the Schwinger functions in the ef- 
fective couplings Ch that is finite at all orders^ see [14] 
for the proof; this is in contrast with the naive perturba- 
tion theory, which is plagued by logarithmic divergencies. 
The boundedness of et makes such expansion meaning- 
ful and it allows one to control the flow of Zh and v^. 
One finds that: (i) lim^_j._oo Vh = 1, i.e., Lorentz invari- 
ance spontaneously emerges; (ii) both Z^ and 1 — v^ 
vanish with two anomalous power laws, see [14]. There- 
fore, if (j)g_j — 0, the dressed propagator has a form sim- 
ilar to Eq.(4), with Z and v replaced by Z(k — k^) and 



w(k- 



-k^J; 



if k is far from the Fermi points, Z(k— k^) and 



w(k — k^) are close to their unperturbed values, namely 
1 and ft. On the contrary, if |k'| < 1, Z(k') - |k'|-'', 

2 

with rj = j^;^ + ■ ■ ■ an anomalous critical exponent 
that is finite at all orders in e^, and w(k') tends to the 
speed of light. Moreover, 1 — u(k') ^ (1 — w)|k'|'', with 
fj — ■^— J- + • • • another anomalous critical exponent. 

The above analysis confirms, at all orders and in the 
presence of a lattice cut-off, the results found long ago 
in [6] , where graphene was described by an effective con- 
tinuum Dirac model with an ultraviolet dimensional reg- 
ularization; the exponents agree at lowest order. Gauge 
invariance implies that the dressed Fermi velocity has 
a universal value (the speed of light) and Lorentz invari- 
ance emerges; this is what happens both in [6] (thanks to 
the use of dimensional regularization) and in the present 
more realistic lattice model, thanks to the exact lattice 
WI Eq.(5). On the other hand, we expect that as soon 
as gauge invariance is broken the limiting Fermi velocity 
is smaller than the speed of light; this is indeed what 
happens in [14], where gauge invariance is broken by the 
momentum cut-off. 

In order to understand which instabilities are likely 
to occur in the system, we compute response functions 
or generalized susceptibilities (which apparently have 
never been systematically computed before, not even in 
the Dirac approximation); more precisely, for 0^^ = 0, 



we compute R^°''{x,j;y,j') 



(") „(") 









with, 



EA^ 



-{0,Sj),a 



'fo^,-M^+s{0fyfi)ds j^ ^^^\ 



(CDW) 

or p'^. 



E.(«^ 



The 



■ - 6+ - b~ 

two latter operators describe, respectively, inter-node and 
intra-node excitonic pairings of K and CDW type; this 
is because the possible presence of a condensate in the 
k = k^ channel for p^ or in the k = channel for 

Px i would signal the emergence of Long Range Order 
(LRO) of K-type (see Fig.l) or of CDW-type (a period- 
2 alternation of excess/deficit of electrons in the sites of 
the A/B lattice). Other relevant bilinears are the Cooper 
pairings, i.e., linear combinations of terms of the form 
at „aZ „ or at- Jj , - , . The large distances asymp- 

X,<T x,-<T x,cr x+(0,(53),cr' ° -' ^ 

totic behavior of the response functions is: 



i?'"'(x,j;0,j) 



g1"'(x) + cos(p+; 



f)G(")(x), 



with |G^" (x)| ^ (const.) |x| ^^ two scaling invariant 
functions (similar formulas are valid for j 7^ j'). In the 
absence of interactions, ^j" — 4, for all a and i; there are 
no preferred instabilities. The presence of the interaction 
with the e.m. field removes the degeneracy in the decay 
exponents: some responses are enhanced and some other 
depressed. It turns out that 



i- 



{CDW) 



= 4-4eV(37r2 



ef^=4-4eV(37r2)+- 



On the contrary, ^2 ^nd ^} are vanishing at sec- 

ond order, while all the Cooper pairs responses decay 
faster than |x|^^. The conclusion is that the responses 
to excitonic pairing of K or CDW type are amplified by 
the e.m. interaction: in this sense, we can say that the 
e.m. interaction induces quasi-LRO of K and CDW type. 
Correspondingly, possible small distortions or inhomo- 
geneities of the K or CDW type are dramatically en- 
hanced by the interactions. For instance, let us choose 
(/)jj as in Eq.(2). The RG analysis can be repeated in 
the presence of a Kekule mass term AoX]s7COs(p^ • 

{x — 6j + Sjg))p^ ■ in the Hamiltonian, which produces 
a new relevant coupling constant, the effective Kekule 
mass. Therefore, the interaction produces an effective 
momentum-dependent gap A(k) that increases with a 
power law with exponent 77^ from the value Ag up to 

A(k±)=Ay(i+''"), ry^ = 2eV(3^^) + --- (6) 

Note that the ratio of the dressed and bare gaps diverges 
as Ao — >■ 0. The enhancement of the dressed gap is re- 
lated to the phenomenon of gap generation in [8], but 
it is found here avoiding any unrealistic large- A^ expan- 
sion. A similar enhancement is found for the gap due to 
a CDW modulation. 

Finally, let us discuss a possible mechanism for the 
spontaneous distortion of the lattice and the opening 
of a gap (Peierls-Kekule instability). We use a varia- 
tional argument, which shows that a Kekule dimeriza- 
tion pattern of the form Eq.(2) is a stationary point of 



the total energy af^ I]^^ >|j + £^o({<^xj}), where the 
first term is the elastic energy and Eo{{(f>gj} is the elec- 
tronic ground state energy in the Born-Oppenheimer ap- 
proximation. The extremality condition for the energy 






where (•) is the ground state av- 
erage in the presence of the distortion pattern {(/'jj}. 

Computing (p^ , ) by RG with the multiscale analysis 
explained above, we find that Eq.(2) is a stationary point 
of the total energy, provided that (fiQ — cqq^ j k + • • • for 
a suitable constant cq and that Aq satisfies the following 
non-BCS gap equation: 



Ao 



9 



dk'- 



z-\k')A{k')\n{k')\^' 



A<|k'|<l 



ki + i;2(k')|r!(fc' +p^)|2 + |A(k')|2 



where A 
-V 



J/ii+v'') 



^0 



and, for A < |k'| < 1, Z(k') - 
Ik'!"", 1 - v{k') - (1 - v)\k'f> and A(k') ~ Aq Ik']-"'". 
In the absence of interactions, Z{k') — 1, v{k') — v and 
the above equation reduces to the free one in [9]. Our 
gap equation is qualitatively equivalent to the simpler 
expression 



1 = 9' 



A^^^i-a- v)pf 



(7) 



from which its main features can be easily inferred. 

At weak e.m. coupling, the integral in the right hand 
side (r.h.s.) is infrared convergent, which implies that 
a non trivial solution is found only for g larger than a 
critical coupling gc] remarkably, gc ~ ^/v, with v the free 
Fermi velocity, even though the effective Fermi velocity 
flows to the speed of light. Therefore, at weak coupling, 
the prediction for gc is qualitatively the same as in the 
free case [9] ; this can be easily checked by noting that the 
denominator in the r.h.s. of Eq.(7) is sensibly different 
from V only if p is exponentially small in v/fj. 

On the other hand, if one trusts our gap equation 
also at strong e.m. coupling and if in such a regime 
ry^ — 77 = Y^p- -I- • • • exceeds 1, then the r.h.s. of Eq. (7) di- 
verges as A — > 0, which guarantees the existence of a non- 
trivial solution for arbitrarily small g; this can be easily 
checked by rewriting Eq.(7), up to smaller corrections, as 



1 



Ai+v-v ^ that is Ao 



.,2(l+,,^)/(r,^-r,-l). 



T}^ — r;— 1 

note the non-BCS form of the gap, similar to the one ap- 
pearing in certain Luttinger superconductors [16]. The 
existence of a non trivial solution for arbitrarily small g 
suggests that strong e.m. interactions between fermions 
enforce the Peierls-Kekule mechanism and facilitate the 
spontaneous distortion of the lattice and the gap gener- 
ation, by lowering the critical phonon coupling gc] this 
is in agreement with the one-dimensional case, where the 
Peierls instability is enhanced by the electronic repulsion, 
see [17]. Note the crucial role in the above discussion 
played by the momentum dependence of the gap term 



(leading to the factor p^~'^'^ in Eq.(7)); on the contrary, 
the growth of the Fermi velocity plays a minor role. A 
similar analysis can be repeated for the gap generated by 
a CDW instability. 

In conclusion, we considered a lattice gauge theory 
model for graphene and we predicted that the electron 
repulsion enhances dramatically, with a non-universal 
power law, the gaps due to the Kekule distortion or to a 
density asymmetry between the two sublattices, as well 
as the responses to the corresponding excitonic pairings. 
Moreover, we derived an exact non-BCS gap equation 
for the Peierls-Kekule instability from which we find ev- 
idence that strong e.m. interactions facilitate the spon- 
taneous distortion of the lattice and the gap generation, 
by lowering the critical phonon coupling. 
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dp 



2^r^ l2;Fp V'k+prA'(k,p)V'k where F^ are matri- 
ces which can be deduced from Eq.(l), with ro(k'^, 0) = 
-il, ri(k^,0) = -aa, r2(k'^,0) = -coai 
Wh is defined in the same way as W , with an extra in- 
frared cutoff suppressing momenta smaller than M^ 

More precisely, ^^[S\^+pKy {'k,-p)S\^]ij = 

9 e^ 



^P aA. M\l 
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